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Abstract 

The vacuum Einstein equations in 5+1 dimensions are shown to admit solutions describing 
naked singnlarity formation in gravitational collapse from nonsingular asymptotically locally flat 
initial data. We present a class of specific examples with spherical extra dimensions. Thanks to 
the Kaluza-Klein dimensional reduction, these spacetimes are constructed by lifting continuously 
self-similar solutions of the 4-dimensional Einstein-scalar field system with a negative exponential 
potential. Their existence provides a new test-bed for weak cosmic censorship in higher-dimensional 
gravity. In addition, we point ont that a similar attempt of embedding Christodoulon’s well- 
known solutions for massless scalar fields fails to capture formation of naked singnlarities in 4+1 
dimensions, due to a diverging Kretschmann scalar in the initial data. 
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I. INTRODUCTION 


In the process of gravitational collapse, singularities can emerge under a broad range of 
circumstances. They are either covered by trapped surfaces (black holes) or being visible to 
far-away observers (naked singularities). Regarding these two outcomes, the weak cosmic 
censorship conjecture [l, 2| claims that generically all singularities should be hidden inside 
black holes. Despite huge effort, the validity of this conjecture still remains one of the most 
important open problems in classical general relativity (see recent reviews 0-^). Neverthe¬ 
less, from a different perspective, signihcant progress has been made on the mathematical 
studies of formation of trapped surfaces in vacuum spacetimes without symmetry assump- 
tion, notably in the recent monumental work by Christodoulou and related development 


15 |. 


In parallel to advancement in 4 dimensions, string theory and brane-world models 


16| 


have greatly promoted the study of general relativity in higher dimensions. By entering 
dimensions higher than 4, novel topologies and horizon geometries can arise, generating a 
large variety of black objects [l7, 18|. Particularly on the dynarnical aspect of the vacuum 


Einstein equations, Choptuik’s type of critical phenomenon 191.120 1 was discovered in collapse 


of 5-dimensional gravitational waves via numerical simulation |2l|]. The spacetime has the 
symmetry of deformed 3-spheres (Bianchi IX) so that Birkhoff’s theorem can be evaded. 
At the threshold of black hole formation, the critical evolution is expected to give rise to 
a naked singularity. Also in 5 dimensions, the Gregory-Laflamme instability 22| of black 
strings has been shown numerically to result in naked singularities through a process of 


self-similar cascading 


23| . Clearly, higher-dimensional spacetimes allow more possibilities 


that cosmic censorship could be challenged. 

In searching of naked singularities, a powerful way to construct solutions in higher di¬ 
mensions is by dimensional reduction j24j. |25| . Originally proposed by Kaluza and Klein, 
the procedure requires that all coordinate dependence of extra dimensions be dropped. The 
spacetime is assumed to be the product of a 4-dimensional Lorentzian manifold and a homo¬ 
geneous space representing extra dimensions. With this approach, the dimensional reduction 
on the vacuum Einstein equations leads to 4-dimensional Einstein gravity coupled with a 
scalar held with or without a negative exponential potential, depending on the geometry of 
the extra dimensions. Under spherical symmetry, the Einstein-scalar held system without a 
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potential (Klein-Gordon) has been thoroughly investigated by Christodoulou. His analysis 


of 


26| established the hrst example of naked singularity 


brmation from regular asymptoti¬ 


cally flat initial data with a suitable matter held (see also ^). This result has been recently 
extended to include the negative exponential potential (Liouville type) by two of us [2^ , and 
new examples that bear similar strength to Christodoulou’s original one were constructed. 

In this paper, we will consider the uplift of these 4-dimensional naked singularity solutions 
and decide whether or not their higher-dimensional descendants can represent formation of 
naked singularity in vacuum gravitational collapse. In particular, we will point out certain 
complications (loss of regularity) that such inheritance can be broken. Our objective is to 
show the following main result: 


In 5+1 dimensional spacetimes with the spatial symmetry of 3“^ x 3“^ orthogonal to a 
radial direction, for the vacuum Einstein equations = 0, there exist nonsingular 
asymptotically locally flat initial data that lead to naked singularity formation. 


The paper is organized as follows. In Sec. II, we present our higher-dimensional metric and 
review properties of the solutions to be embedded. For the negative exponential potential, 
we examine the 6-dimensional solutions in Sec. Ill from the aspects of apparent horizons, the 
Kretschmann scalar, asymptotic behavior, and homothetic symmetry. Sec. IV comments on 
a similar treatment of Christodoulou’s solutions in 5 dimensions. Concluding remarks are 
made in Sec. V. We use the same notations of 2^ throughout the paper. For the conceptual 


framework involved in this work, one may refer to 


29l | (Chs. 1 & 2), which contains a quick 


survey of general relativity from the point of view of the Cauchy problem. 


II. THE LIFTED SOLUTIONS 

We start with the following warped product metric (hatted) in 4 -|- n dimensions with a 
parameter k\ 

ds^ = —g{x)g{x) — 2g{x) du dr -|- r'^dVfl, (2) 

- r 

(t){u,r) = h{x) — Khi{—u), K > 0, x = — >0, m < 0. (3) 
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Here the 4-dimensional metric ds^ (unhatted) is written in retarded Bondi coordinates, and 
dfi^ denotes the metric of a unit 2 -sphere and ds^ a compactifying space fC having a constant 
Ricci scalar Rn = —4Vo with Vq = 0 or —1. For simplicity, we assume /C = S'” or T”. By 


virtue of consistent Kaluza-Klein dimensional reduction on ds^ 


24 


25 


30l-l32| , one can verify 


that the vacuum Einstein equations for ([T]) correspond to the equations of motion from an 
effective action on the metric ds^ and the scalar held 0 (dilaton): 


£ = (f - - V{4 >)), 


V = V'oexpl —) , 


K j 


K = 


> 0 , 


if Ro = 0 {Rn = 0, /C = T”), 


V 


n 

n+2 


(4) 

(5) 


< 1, if Vo = -1 (Rn = 4, /C = R 

Hence the effect of the extra dimensions is tantamount to actuating a scalar held in 4- 
dimensional Einstein gravity, of which the potential V{(j)) depends on the choice of /C. In¬ 
troducing new unknown variables ( ' = d/dx) 

C(s) 




X 


7 (s) = xh', X = e^ > 0 , 


( 6 ) 


2(1-cw) s ' 

we can reduce the Einstein equations to a 3-dimensional autonomous system of hrst-order 
ordinary diherential equations (' = d/ds = xd/dx) 

2/(7+ k)(37C-27 + 0 


y = 


i-C 


C = - [2(7 + k)0 1 - + [(7 + nf + 1- c, 

7 = 7^ + ( k “0 2 k ) 7^ + ^ ^ 


1 - 2C 

Additionally if Vq 7 ^ 0 , the function h is subject to 


k(1-2C) 


1 - - 
y 


(7) 

( 8 ) 
(9) 


(7 + 


2C .^2 

= 1 + K 


1 - 2Vox^e‘^~^/^ 


1-C y 

It is straightforward to verify that flTl fTOj) together imply the equation 7 = xh' in 


( 10 ) 


If 


Vo = 0 as in the case of massless scalar helds (/C = T”), the equation flTOl) without h then 
serves as a constraint to the system (1319]), generating a 2 -dimensional submanifold in the 
phase space. 

With a combined analytical and numerical approach, the dynamical system (I3l9|) and the 
associated 4-dimensional metric ([2|) have been systematically studied in 2^, which extends 


the earlier results by Christodoulou 


26l | and Brady j27| on the special case of massless scalar 
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fields. It is established in 4 dimensions that solntions describing naked singnlarity formation 
(see Fig. [T]and[2]) are represented by the integral curves of ([319]) that run continuously from 
the initial starting point 


0-. y = I, C = 0 , 7 = 0 , s =-cx), 


( 11 ) 


to the ending points (forming a line) 


'Poivo) ■■ y = Vo 


1 


if Vi) = 0, 
if Vi, = -1, 


c = 1, 7 = -K, s = +CX). (12) 


The condition (HD is imposed by regularity at the symmetry center r = 0, u < 0 (x = 0). 
Requiring the solution curves to terminate at the line Vo{y) ensures that the future light 
cone M = 0 of the curvature singularity can escape from the center. Lack of an apparent 
horizon in the spacetime is guaranteed by 


= y{s) ^ 0, -oo < s < +oo. 


(13) 


Fulhlling all these conditions, a 2-parameter family of naked singularity solutions can exist 
provided 

0 < < 1/3 for Vo = 0, (14) 


or 


0 < < 1 for Vi) = —1. 


(15) 


Note that among these two possibilities, k is a free parameter in the former but hxed by the 
potential (0]) in the latter (cf. ([5])). Furthermore, for each solution in the family, there lies 
a midpoint s = s* such that the integral curve reaches a singular point of the ODEs: 


Ps{y*) ■ y = y* 




> 


l + K^ 


if Vi, = 0, 
if Vi, = -1, 


c = 


2 ’ 


7 = 


^-y* 

y*tP 


s — s*. 


(16) 


It bears the signihcance that the past light cone of the central singularity is given by —r/u = 
X* = exp(s*). By identifying the asymptotic behavior of the solutions as s —)■ 


deduce that in the limit u —)■ 0, a curvature singularity forms at the center 
Ricci scalar obeying 


26 


oo, one can 


28 | with the 


R — 


2 ( 2/0 - 1 ) ^ - 2 - AyoK 




, as M —)■ 0 (x —)■ -|-oo). 


(17) 
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Meanwhile, the scalar held diverges logarithmically at r = 0: 


-)■ 


—Klnr if Vo = 0, 

-Klnr + f if Id) =-1, 


as M —>■ 0. 


(18) 


In fnture sections, we will determine whether this basic pictnre of 4-dimensional naked 
singularity spacetimes as illustrated in Fig. [T]may carry over to higher dimensions. 



FIG. 1. A spacetime diagram of the naked singularity formation near the symmetry center. The 
self-similar spacetime is obtained by solving the autonomous system ([7][9]) from x = 0 to x = -|-oo 
(the coordinate line x = —rju sweeping from r = 0 to u = 0). The line x = x* (x = -|-oo) marks 
the past (future) light cones of the central curvature singularity at r = 0, « = 0. There is no 
apparent horizon. The curvature invariants (the Ricci and Kretschmann scalars) are well-behaved 
before n = 0 at the center and beyond. 


III. NAKED SINGULARITY FORMATION IN 6 DIMENSIONS 

Returning to the metric ([T]), we hrst consider the uplift of the 4-dimensional naked sin¬ 
gularity solutions with Vq = —1. We will show that the resulting spacetimes inherit major 
properties of the lower-dimensional ones, and thus can represent formation of naked singu¬ 
larities as well (cf. Fig. [1]). To make our calculation concrete, we focus on the simplest case 
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FIG. 2. Numerical integration of the autonomous system (i7ra with Vq = —1 and k = l/\/2. 
Integral curves from the initial point O {x = 0) first reach the singular curve Vs{y) (x = x*, 
( = 1/2), and then they bifurcate and continue to move towards the line Voiy) {x = +oo). The 
plot shows such bifurcation drawn from one sampling point on Vs{y)- 

with /C = K = l/\/2 and take 

dsl=2 = (19) 

with the metric of a unit 2-sphere. An extra factor 1/2 appears above so that the Ricci 

scalar of JC follows Rn =2 = 4 (cf. (|5])). Under the assumptions, the metric ([T]) acquires the 
spatial symmetry of x 


A. Apparent horizon 


The radial null geodesics of the metrics ([T]) and 
outgoing null rays and 


dr 

du 


I 

2’ 


are both determined by n = const for 


( 20 ) 


for incoming null rays (for construction of double-null coordinates, see 2^, Sec. 4). They 
correspond respectively to the future-directed null vector helds I and n: 

2 


1 9 ^ 

I = \ -Or, 

V 9 


n = 




du 


'9 Q 

Ut') 

9 


ds^(/, n) = —2, 


( 21 ) 
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which are normal to a 4-dimensional spacelike surface S = x parameterized by r, m = 
const. We introduce the null second fundamental forms y and y of the surface S, which are 
defined by 

x(X, y) := ds^iVxn, Y), x(X, F) := ds^Vxl, Y), (22) 

for two arbitrary vectors X and Y tangent to S. Their traces measure the expansions of 
the two radial null geodesic congruences along I and n. In this setting, an apparent horizon 
(marginally outer trapped surface |33M37l|) consists of a collection of E’s such that 


trx|E = 0. trx|j<0. 


(23) 


Thus it boils down to calculating these two quantities. The result reads 


ry g 


trx = . = —W-, 

r \l g 


ig 


(24) 


for general input of = g{u,r), g = g(u,r), and 0 = (j){u,r) to the metric ([T]) with (IT^ . 
Remarkably, the same expressions also hold for the metric ds^ (E = S^) itself, as well as 
when JC = (E = x T^, cf. (IT7|) i. This is due to certain key cancelations enabled by the 
Kaluza-Klein reduction. Hence identically to the situation in 4 dimensions (cf. flT^ b absence 
of the apparent horizon is signaled hj g/g = y ^ 0 in the lifted spacetime. For similar 
discussions on such invariance of trapping within the Kaluza-Klein dimensional reduction. 


one can see 


38 


B. Kretschmann scalar 

Much like its 4-dimensional counterpart, the lifted spacetime can maintain nonsingular 
for M < 0, both at the center and beyond. A scalar-valued curvature singularity first occurs 
at r = 0, M = 0. To see this, one should examine the Kretschmann scalar K = Rfiupa 
particularly at three key locations in the spacetime, i.e., a; = 0, -|-cx) and x^. In terms of ?/, 
Q and 7 , the full expression of K is given by 

K = 12{ [ - 6(2C - 1) V - 2\/2(9C - 4)(2C - 1 ) 7 ' - C(25C - 12 ) 7 ' 

-4^/2C^ - 9C' + 16C - 8]|/2 + [ - 8(2C - 1)(C - 1)7^ 

-4^/2C(C - 1)7 + 16(C - i?]y - 8(C - 1)'} 

/r\C - 1) [ 27(7 + V^)(2C - l)y - (C - 2 ) 2 / + 2(C - 1)], (25) 
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where we have used the equations ([319]) to eliminate all the derivatives of y, ( and 7 . 

Near the initial point O with a; = 0, the system ([319]) admits the following Taylor series 
solution 


y = 1 — ctx^ — 2ctx^ + 

( = 2x — 4a;^ + (2 ct + 8 — k‘^)x^ + 0{x‘^), 


(26) 


-y = KX + [ — + k]x'^ + 


3ct , ,, 

— + {ct + 1)k-- 


X^ + 0 ( 3 ;“^), 


with a parameter ct < k^/3 (ct = k^/ 3 for Vq = 0). Plugging it into (l25D . we obtain 


K = 


18(1—4cr)^ ^ 2\\ 1 

-j—(1 + 2i + 0(1 )). CT<Y = Q- 


(27) 


which stays hnite and continuous at r = 0 
Voiyo) with x = +cx), the solution satishes 

1 2 


or hnite u < 0. Similarly, near an ending point 


28| 


y^yo> T—^ (C - i)e^^ ^ C 2 < 0 , 

L Kj O 


(7 + , 

V yoi^ J 

where C 2 is a constant. These asymptotic limits result in 


as s —)■ +CX), 


K 


122/0 


as s ^ +CX) (u —>■ 0—). 


( 28 ) 


(29) 


{3vo-2)r^’ 

Hence, the scalar K diverges at r = 0, but remains bounded elsewhere along u = 0 (the 
Cauchy horizon). Related to this, the same fall-off behavior also holds for other hxed u < 0: 


K = 


122 / 0 ^ 


-|- o(r ) —)■ 0 , as r —)■ -|-cx). 


( 30 ) 


(St/o - 2)r-‘ 

Now move on to x = x*. To capture the asymptotic behavior of the solutions, we 
introduce a new independent variable t to the system ([319]) such that 


ds 

= - 2 ( 1 - 20 . 


(31) 


n terms of f, the singular point Vs{y*) turns into a critical point of the transformed equations 
281]. Therefore, we can perform standard linearization with 


1 1 ~ 2 /* 1 

y = y^+Xi{t), ( = -+x 2 {t), 7 = ,...3 +a(3(C, y*> 


2 y^K 

The eigenvalues of the linearized system are 


1 + 


(32) 


Ai = A 2 = 1 — 


A.s — 0 . 


( 33 ) 
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The zero eigenvalue reflects the fact that Vs{y) forms a curve. Given that n = l/\/2, the 
first two eigenvalues coincide: 


Ai = 


A 2 — 


1 

2‘ 


Accordingly, the solution in a neighbourhood of s = s* {t 
Sec. 2) 


—cxd) are of the form (cf. 


(34) 


m 


y 

c 

7 


4(5?/* -4)(3?/* -4)ai ^2 , 

y* -e / H-, 

y* 

^ + aie*/2 + • • • , 


2^/2 {I-y,) 


+ 


4\/2(3?/* - 4)ai 
0,2 - 1 

y* 




(35) 

(36) 

(37) 


where Oi 2 are two parameters and higher-order terms in (ri^ 2,3 are omitted. Consequently, 
we obtain 

12y2 

^ ^ (3?/* -2)r2 ’ - 00 ), (38) 

which resembles the limit fl2^ . Further numerical tests by us agree with fl38|l . 

For other fixed values of x > 0, the equation fl2^ implies that 


iF cx — —)■ 0, as r —)■ -f cx). 


(39) 


where the boundedness of the proportionality coefficient can be verified by numerical calcu¬ 
lation for individual solutions. 


C. Asymptotic local flatness 

With K = l/\/2, the limits in (l28ll imply that 

99 4c| g 2 c 2 \/3?/o - 2 

X ?/o I/O -2 m 

Introducing new coordinates 

= r, i — r = u = 2 c 2 \/—u < 0, 


as X —)■ -|-cx). 


(40) 


(41) 
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( 42 ) 


thus for finite u < 0 with r (equivalently x) being large, we have 


= e^'^ (^—g{x)g{x) du^ — 2g{x) dudr + + 

2^/3yo^^ f c^dn^ 2c2dMdr\ ^ 

J 




■ dn^ 


2/0 

2\/dyo — 2 


+ + 


—u 


2/0 

2\/dyo — 2 


(du^ + 2d-udr) + dhl^ + 


1 


2/0 


{-dP + dh^) + ^ 


Avo 

1 


\/32/o - 2 


- 2 
d02 


- 2 
df2^ 


dff' 


(43) 

(44) 

(45) 
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Then note that by stereographic projection, r^dhl^ = (dx^ + d^^)/[l + (x^ + y^)/4r 
dSP + d^^ for local x and y as r ^ +cx3. Combining these approximations with the decaying 
Kretschmann scalar in fl30l) . we comment that the lifted spacetime is asymptotically locally 
flat with non-trivial topology x at spatial inhnity. 


D. Homothetic Killing vector 

With the assumption of flT^ and k = l/A, one can check that the lifted metric ds^ 
{.9A possesses the same homothetic Killing vector ^ as ds^ {.9A l26|-l28|: 

udfj^ T 'fdri L^9iiu ^^9(ii> 9jj.i'- (^6) 

Therefore, the 6-dimensional naked singularity spacetimes with spherical extra dimensions 
are continuously self-similar. 

IV. COMMENT ON 5 DIMENSIONS 


From the previous discussion, one would anticipate that lifting the solutions with Vq = 
0 for massless scalar helds might also generate new naked singularity solutions in higher 
dimensions (see ^ for discussions on lifting solutions with k = 0). However, we will show 


that this turns out to be false. The reason lies in that, unlike the 6-dimensional example, 
the Kretschmann scalar blows up at x = x*, i.e., along the past light cone of the central 
singularity. Hence in such a spacetime, a curvature singularity can pre-exist before the 
central singularity emerges later at u = 0. 

Again for concreteness of calculation, we focus on the simplest case /C = T^(= S*^) in 5 
dimensions and take 

d4=i = (dx^)2 (47) 
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with the local coordinate on /C. Using (fTOj) with Vq = 0, we can remove the unknown y 
in the equation ([9]) and obtain a reduced system for ( and 7 . It can also be verihed that 
the equations flSl lTUD together imply ([7]). To solve this 2-dimensional autonomous system 


in a neighborhood o 


and obtain 


the singular point one can apply the same procedure as in 


3 . 


Sec. 2 ; for a comparison of notations, see the Appendix of 28|) 


C = ^ + + • • • , 

1 

7 = - + 

K 


kH , 4a;(l + K^) 

■ ^k3(1-2k2) + ’ 


( 48 ) 

( 49 ) 


where the variable t is dehned by flHT]) and 01^2 are two parameters. Owing to these asymp¬ 
totic solutions for t —00 {s —)■ s*), the scalar K (written in terms of C, and 7 without 
derivatives) diverges according to 


lim 


(X 


t^ — OO 


^4-4k/x/3’ 


0 < K < 

\/3 


( 50 ) 


One can also conhrm the divergence of K at x = via numerical integrations along indi¬ 
vidual solution curves. 

As an aside, we note that the lifted metric ds^ with f|T71) can no longer admit the ho- 
mothetic Killing vector ^ = udu + rdr unless n = 1/ \/3, of which the value lies outside the 
required bound 0 < < 1/3 for naked singularities in 4 dimensions. 

At first glance, the inconsistency of spacetime regularity (cf. (ISOD i in the Kaluza-Klein 
dimensional reduction may appear puzzling. One physical way to resolve it is by changing the 
conformal frame for ds^ 41|, which means that one should instead look at the 4-dimensional 


projection of the full metric ds^, i.e., exp(20/\/3)ds^. For this re-scaled metric with the 
dilaton (j) involved, direct calculation shows that the Kretschmann scalar diverges at a; = x* 
{K continuous for ds^ except at r = 0, u = 0; iF cx 1/r^) in accordance with its 5-dimensional 
counterpart. 


V. CONCLUDING REMARKS 

In the context of the Kaluza-Klein dimensional reduction, we have constructed a 2- 
parameter family of 6 -dimensional naked singularity solutions by lifting the 4-dimensional 
ones in the Einstein-scalar held system with a negative exponential potential. Despite its 
success, the same construction does not come through in the case of a vanishing potential 
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due to a diverging Kretschmann invariant along the incoming null ray x = (hence in 
the initial data), which raises a caveat for utilizing such lifting/embedding. Regarding 
the 6-dimensional solutions, it may draw some concern that the spatial inhnity has a non¬ 
trivial topology 5^ X 5^. However, we should emphasize that the local formation of a 
singularity near the center does not depend on the spacetime region far from it. Hence it 
would be interesting to investigate if the central region of the spacetime can be truncated 
and matched to an appropriate surrounding region that is, for instance, asymptotically 
Minkowski (S^) or asymptotically Kaluza-Klein (compact extra dimensions). Related to 


this, we recall that such truncation is necessary for Christodoulou’s solutions 


26| . in order 


to obtain asymptotically flat initial data. As a hnal remark, it should be noted that though 
we have restricted ourselves to 6 and 5 dimensions, our treatment can be generalized to 
other values of 4-|-n > 6 (certain complications may arise). The results with detailed proofs 
will be reported elsewhere. 
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